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SECTION  I 
Introduction 

In  Ref  [l]  and  [2]  the  stability  of  the  lower  branch  solution  of  the 
Falkner-Skan  similar  boundary  layer  equations  was  investigated.  These 
velocity  profiles  possess  one  inflection  point  and  give  rise  to  the  "Rayleigh 
Instability".  The  analysis  of  this  instability  proved  extremely  useful  in 
interpreting  self  excited  oscillation  occuring  in  cavities,  over  spike  tipped 
bodies  and  in  inlets  (Ref.  [3,  4,  5,  6,  7,  8,  9]). 

Other  classes  of  self-excited  oscillations  have  been  observed  in  jets 
(e.g.  edge  tones)  and  in  the  wakes  of  bluff  bodies  (e.g.  periodic  shedding 
of  vortices  behind  cylinders).  The  velocity  profile  for  this  class  of  flows 
possess  two  inflection  points  which  give  rise  to  two  different  modes  of 
instability  (Ref.  [lO,  1 l] ) .  To  assist  in  the  interpretation  of  these 
observed  instabilites  it  was  felt  useful  to  further  investigate  the  stability 
features  of  compressible  wake  and  jet  profiles.  For  this  reason  eigenvalue 
solutions  for  a  series  of  typical  profiles  were  computed  for  the  following 
types, 

2 

(a)  Symmetric  jet  U  =  sech  y 

2 

(b) .  Symmetric  'wake  U  =  sech  y 

3  2 

(c)  Anit-symmctric  U  =  —  sech  y  tanhy 

(Combined  wake  and  jet) 

(d)  Asymmetric  jet  0  -®<y<-2.5 

.23529(y+2.5)2  -2.5<y<-.8 

U  =  1-  .5  y2  -.8<y<1.25 

1 .7857 (y— 1 . 6)2  1.25  <y<1.6 

0  1.6  <y<® 

(e)  Asymmetric  wake  U  ■  -U  (of  case  d) 

The  results  of  the  stability  analysis  are  compiled  and  catologcd  to 
permit  our  conclusions  regarding  the  behavior  of  these  flows. 


Objective 

The  objective  was  to  determine  the  amplification  factor,  disturbance 
propagation  speed  and  wave  number  for  typical  velocity  profiles  with  two 
or  three  inflection  points  at  various  Mach  numbers.  It  was  anticipated  that 
some  overall  characteristics  for  wake/jet  flows  could  be  deduced  from  these 


series  of  calculations. 


SECTION  II 


Governing  Equations 


In  this  report  we  study  the  stability  of  compressible  wakes  and  jets 
in  two  dimensional  flows.  Let  u  represent  the  velocity  component  in  the 
x  direction  and  v  the  velocity  component  in  the  y  direction,  p,  p  and  T 
are  pressure,  density  and  temperature  respectively. 

The  basic  oquaitons  are 


(2.1) 


3u  .  3u  ,  3u  1  9p 

- —  +  u  — —  +  v  - —  =  -  —  — c- 

3t  3x  3y  p  3x 


(2.2) 


3v  .  3v  ,  3v  .  1  3p 

—  +  u  —  +  v  — —  - - 

3t  3x  3y  p  3y 


(2.3) 


1 £ 

3t 


+ 


ul£+vlEei£[iP  + 

3x  3y  p  L  3t 


3p 


9p 


U  3^+  ] 


(2.4) 


Eliminating  p  between  equations  (2.1)  and  (2.4)  we  obtain 
3p  .  3p  .  3p  .  ,3u  .  3v. 

at  +  u  al +  v  a?  +  YP<S  +  a7>  "  0 

Equations  (2.2),  (2.3)  and  (2.5)  have  a  steady  state  solution 
u  **  u(y),  v  ••  0,  p  =  p  =  constant 
we  assume  the  time  depndent  perturbed  flow  as  [l2,13] 


u  =  u(y)  +  u' (x,y,t) 
v  -  v' (x,y, t) 
p  -  p  +  p'(x,y,t) 


(2.5) 

(2.6) 

(2.7) 

(2.8) 
(2.9) 


Substituting  these  values  of  u,  v  and  p  in  equations  (2.2),  (2.3)  and  (2.5); 
and  retaining  only  linear  terms  in  u' ,  v’  and  y>'  we  obtain 


I 


|h'  +  :|H'  +  v’  |5  0 

3t  3x  3y  -  3x 


3v*  -  3v  1  3 £'  _ 

—  +  u  —  +  —  — =  0 
3t  3x  -  3y 


3p'  ,  -  3p'  -  /3u '  ,  3v\ 

ST  +  “  to  +  ™  (3iT  +  3y  >  '  0 


We  seek  the  periodic  solutions  of  the  form 
ia(x-ct) 


(2.10) 


(2.11) 


(2.12) 


o'  =  u(y)  e 


v'  =  v(y)  e 


ia  (x-ct) 
ia  (x-ct) 


(2.13) 

(2.14) 

p«  =  p(y)  (2.15) 

*  A  * 

where  u,  v,  p  are  complex,  c  is  a  complex  constant  and  a  is  a  real  constant. 

Substituting  (2.13),  (2.14)  and  (2.15)  into  equations  (2.10),  (2.11) 
and  (2.12)  we  obtain 

(2.16) 


ia(u  -  c)  u  +  u  v  =  -ia  ^ 
y 


ia(u  -  c)  v 


1  ' 
T  P, 


(2.17) 


ia(u  -  c)  p  *=  -  YP  (iau  +  v  )  (2.18) 

We. eliminate  p  and  u  from  the  above  equations,  use  the  relation 
p  ■*  p  R  T  and  obtain 


(u  -  c) 

V 

-  U  V 

y 

y 

,  2 

Y  R  T  - 

(Z 

-  c) 

q^(u  -  c)  u 


Y  R  T 


Now  using 

Y  R  T 


1+.2H2  (1  -  Z2) 


(2.19) 


(2.20) 


and  doing  some  calculations  we  obtain 
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S' 


(u  -  c)  Vy  -  »y  V 
(1+.2M2)  -  M2  (u  -  c)2  y 


=  a2  (1  +  . 2M2)  (»  -  c)  v  (2.21) 


where 


o 


M 


2 


1+.2  M2 
00 


(2.22) 


and 


-2 

u 


1+.2  M2(1-u2) 

CO 


(2.23) 


If  we  write 

2-12-  A  - 
(1  +  .2M  )  -  M  (u  -  c)  =  g,  and  replace  v  by  )  and  u  by  1)  in  equation 

(2.21)  we  obtain 


(U  -  c)  *y  -  Uy  $ 
fi 


a2(l 


,2M2) 


(U  -  c)4> 


(2.24) 


For  boundary  conditions  we  assume  that  for  unbounded  flows  the  initial 

disturbances  die  down  at  far  from  the  disturbances.  Therefore  we  get 

*  (-“)  =  0,  *(“)  *  0  (2.25) 

For  fined  wave  numbers  (a  =  constant)  equations  (2.24)  and  (2.25)  is 

an  eigenvalue  problem.  |  is  eigenfunction  and  c  is  eigenvalue. 

We  solve  this  eigenvalue  problem  for  the  following  velocity  profiles 
2 

U(y)  =  sech  y,  symmetric  jet  (2.26) 

2 

U(y)  *  -sech  y,  symmetric  wake  (2.27) 

3  2 

U(y)  *  j  r'j  sech  y  tanhy,  anti-symmetric 

(combined  wake  and  jet)  (2.28) 


u(y) 


w  2 
•23529(y+2,5) 

1  -  .5y2  2 

1.7857(y-1.6) 

0 


— °°  <  y  <  -2.5 

-  2.5  <  y  <  -  .8 

-  .8  <  y  <  1.25 
1.25  <  y  <  1.6 
1.5  <  y  <  » 


asymmetric  jet 


(2.29) 


Asymmetric  wake 


SECTION  III 


Numerical  Procedure 


Eigenvalues  of  4>  were  determined  by  a  shooting  mfcthod  [1]:  starting 

with  boundary  conditions  at  ymin»  integrating  over  the  range  of  y,  and 

comparing  the  result  with  the  outer  boundary  condition,  namely  $  =  0  at 

Y  .  The  process  involved  minimization  of  the  error  caused  by  the 
max  r 

2  2  2 

deviation.  This  was  chosen  to  be  the  square  of  the  norm  of  4>,  |$|  =  4>  +4^. 
The  integration  was  done  using  a  fourth-order  Runge-Kutta  method. 

Boundary  conditions  at  ymin  were  determined  by  observing  the  behavior 


of  (2.24)  as  y  -*■-  <».  The  equation  reduces  to 
4  =  a  4 

yy 

Since  we  desire  4>  (— °»)  =  0,  we  choose 

l«|y. 


4>(y  .  )  -  e 
mm 


min 


*'<W  ■  bio 


bkin 


(3.1) 


(3.2) 


as  our  boundary  conditions. 

The  method  of  finding  eigenvalues  utilized  the  same  minimization 
routine  as  in  previous  investigations  [l,2].  The  user  provides  a  storing 
guess,  for  c  in  the  case,  and  the  routine  begins  by  searching  along  a  con¬ 
stant  line  of  c^  with  increasing  steps  until  the  error  begins  to  increase. 

It  then  uses  the  last  three  calculated  values  to  determine  a  parabola,  with 
the  c^.  value  at  the  vertex  used  as  the  new  approximation.  Then  this  value  of 
cr  is  held  constant  and  a  search  along  a  line  of  changing  c^  is  carried  out. 
After  a  new  relative  minimum  is  found,  the  quadratic  approximation  is  used 
to  determine  a  new  value  for  c^.  The  third  step  involves  searching  the  line 
connecting  the  original  guess  and  the  new  point  in  the  same  manner.  If  the 
error  is  not  less  than  a  preset  limit,  here  10  **,  the  routine  starts  again 
with  the  latest  value  used  in  place  of  the  original  guess. 


SECTION  IV 


RESULTS 


The  eigen  value  problem  represented  by  (2.24)  and  (2.25)  was  solved 
numerically  for  the  velocity  profiles  given  by  (2.26),  (2.27),  (2.28),  (2.29) 
and  (2.30).  The  results  are  tabulated  for  a  wide  range  of  wave  numbers  (a) 
and  Mach  numbers  (M^) .  The  instability  characteristics  for  a  ay  metric  jet, 
asymmetric  jet  and  anti-symmetric  jet  are  given  in  tables  (la)  (3h).  For 
=  0  there  values  agree  with  those  given  in  [10  and  ll]. 

The  velocity  profiles  are  plotted  in  Figures  (1-3).  Tt  ies  of 

a,  versus  c^,  and  a  versus  and  c^  versus  c^  are  plotted  in  Figures  (4  -  19). 
<|>,  u  and  p  are  plotted  for  some  special  values  of  a  and  c  in  Figures  20  a  to 

21  c,  and  mangitudes  and  phases  of  <J>,  u  and  p  are  plotted  in  Figures  22a  to 

22  c,  and  28  a  to  28  c. 

Solutions  were  obtained  with  convergence  error  criteria  of  at  least 
10‘6  for  all  cases. 


SECTION  V 


Summary 

The  stability  of  compressible  inviscid  jets  and  wakes  has  been 

v 

investigated  by  utilizing  the  linearized  equations  resulting  from  a  small 
perturbation  analysis.  The  resulting  eigenvalue  prohlems  were  solved 
numerically  for  various  wave  numbers  (2)  and  Mach  numbers  (M  )  for  differ- 
ent  velocity  profiles.  In  the  cases  of  symmetric  jets  and  wakes  and  that 
of  asymmetric  jets  and  wakes  we  found  two  propagation  modes  corresponding 
to  two  inflection  points.  The  sinuous  mode  for  even  eigenfunctions  and 
varicose  mode  for  odd  eigenfunctions. 

In  varicose  modes  the  magnitude  of  amplification  decreased  as  Mach 
number  (M  )  increased  and  the  flow  became  completely  stable  at  =  2. 

In  sinuous  modes  the  amplification  did  decrease  a  little  with  the  increase  of 
Mach  number  but  we  did  not  find  any  upper  limit  in  Mach  number  above  which 
the  flow  was  completely  stable. 

In  the  case  of  anti-symmetric  profile  there  are  three  modes  correspond¬ 
ing  to  the  three  inflection  points.  Two  propagating  modes,  one  propagating  to 
the  right  and  the  other  propagating  to  the  left;  and  one  standing  mode.  The 
magnitude  of  amplification  for  propagating  modes  decreased  as  the  Mach  number 
increased,  and  completely  died  down  at  Mach  number  of  1.5.  On  the  other  hand, 
we  could  not  find  an  upper  limit  of  Mach  number  for  the  standing  mode  above 
which  the  flow  was  completely  stable.  The  authors  believe  that  these  results 
will  be  useful  for  analyzing  aerodynamic  instabilities  encountered  in  wakes 


and  jets. 
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INSTABILITY  ClIAKACTEKiSTICS  KOU  TDK  SYMMKTM1C  JET 

TABLE  1. 

la)  M^.0.0,  SINUOUS  MODE 


.061256 

.137549 

.207237 

.266554 

.316088 

.357248 

.392290 

.422860 

.450120 

.474924 

.497882 

.519408 

.539851 

.559444 

.578370 

.596770 

.614738 

.632352 

.649655 

.666667 


=scch2y 


c. 

l 

.119380 

.205118 

.241188 

.249623 

.244302 

.231763 

.215421 

.197142 

.177992 

.158612 

.139405 

.120635 

.102479 
*■  .085063 
.065476 

.052789 

.038056 

.024322 

.011026 

.282007(10) 


t-  00  0> 


lb)  =1.0,  SINUOUS  MODI'. 


c 

c . 

a 

r 

l 

.1 

.067994 

.128268 

.2 

.156481 

.216384 

.3 

.233560 

.245433 

.4 

.296861 

.246748 

.5 

.346972 

.236112 

.6 

.387995 

.219914 

.7 

.422773 

.201173 

.8 

.453158 

.181434 

.9 

.480444 

.161569 

1.0  .505465  .142098 

1.1  .528781  .123328 

1.2  .550770  .105484 


1.3 

1.4 

1.5 


571686 

591593 

610908 


088689 

073020 

058512 


1.6  .629391  .045172 
-1.7  .647185  .032987 
i.«  .664327  .021922 


1.9 

2.0 


680823 

696684 


011935 

022972 


lc)  M  -2.0,  SINUOUS  MODE 

CO 

c 

c. 

a 

r 

l 

.1 

.086251 

.149243 

.2 

.203675 

.235800 

.3 

.293797 

.247300 

.4 

.359668 

.234900 

.5 

.410110 

.214300 

.6 

.451142 

.191700 

.7 

.486228 

.169100 

.8 

.517361 

.147500 

.9 

.545685 

.127500 

1.0 

.571836 

.109323 

1.1 

.596155 

.093041 

1.2 

.618823 

.078622 

1.3 

.639958 

.065944 

1.4 

.659655 

.054849 

1.5 

.678011 

.045186 

1.6 

.695111 

.036856 

1.7 

.710884 

.029902 

1.8 

.724660 

.023967 

1.9 

.737094 

.017381 

2.0 

.749815 

.011191 

2.1 

.761940 

.006051 

(OOtxJ  9  Ut  fr  U  N 


a 

.1 


1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 


M  =3.0,  SINUOUS  modi; 

00  1 

c 

r 

c. 

i 

. 112176 

.172854 

.259379 

.247000 

.355744 

.237003 

.422588 

.214158 

.474158 

.186233 

.517294 

.159762 

.555257 

.136500 

.589379 

.117193 

.619734 

.102317 

.645257 

.091204 

.665674 

.081286 

.683340 

.070539 

.700638 

.059729 

.717718 

.050427 

.733757 

.042849 

.748628 

.036578 

.762489 

.031398 

.775722 

.027867 

.785647 

.027098 

.791696 

.022900 

le) 


1^4.0,  SINUOUS  MODE 

c 

c. 

r 

X 

.143092 

.193015 

.311005 

.249100 

.409095 

.225500 

.477423 

.191509 

.532312 

.159900 

.579928 

.135133 

.620580 

.119050 

.650447 

.109428 

.671343 

.098700 

.689858 

.085042 

.709335 

.071677 

.728751 

.061123 

.746833 

.053174 

.763496 

.047336 

.778596 

.044463 

.788371 

.043617 

.795081 

.038781 

.804006 

* 

.032273 

813824 

823122 


027343 

023604 


If)  Moo=0 . 0 ,  VAIilCOSE  MODE 
c 

C. 

a 

r 

l 

.05 

.862061 

.030296 

.10 

.867554 

.108700 

.20 

.796327 

.121800 

.30 

.759672 

.114815 

.40 

'.733435 

.102812 

.50 

.713113 

.088200 

.60 

.697180 

.071556 

.70 

.684964 

.053968 

.80 

.676071 

.035862 

.90 

.670036 

.017700 

1.00 

.666667 

.000000 

lg)  Hb- 1.0,  VARICOSE  MODE 
c 

c  ■ 

a 

r 

1 

050 

.814895 

.011590 

100 

.834826 

.151395 

200 

.765617 

.075559 

300 

.737390 

.056628 

400 

.721227 

.039567 

500 

.710643 

.022751 

525 

.708624 

.018530 

550 

.706811 

.014300 

575 

.705210 

.010064 

600 

.703712 

.005774 

625 

.702614 

.001584 

650 

.701492 

.000049 

lh)  M^-2.0,  vaiucose  mode 

c 

c. 

a 

r 

l 

54 

.988650 

.000784 

55 

.980628 

.001722 

56 

.973019 

.002772 

57 

.965803 

.003849 

58 

.958940 

.004919 

59 

.952416 

.005956 

60 

.946200 

.006928 

61 

.940272 

.007830 

62 

.934607 

.008647 

63 

.929187 

.009374 

64 

.923988 

.010011 

65 

.918995 

.010557 

66 

67 

68 


914191 

909559 

905084 


.011011 
.011378 
■ .011660 


INSTABILITY  CI1ARACTER1 ST ICS  FOR 

r\ 

THE  ANTISYMMETRIC  JET  U  *=  'j  ft  scch  ytanhy 

TA^Le  2. 

2a)  M00=0 . 0  ,  PROPAGATING  MODE 


a 

c 

r 

c . 

l 

.05 

.920689 

.090656 

.10 

.879349 

.116531 

.20 

.820211 

.136200 

.30 

.771224 

.138417 

.40 

.730500 

.130034 

.50 

.699088 

.115459 

.60 

.675954 

.098057 

.70 

.659475 

.079866 

.80 

.648131 

.061943 

.90 

.640730 

.044786 

1.00 

.636391 

.028606 

1.10 

.634469 

.013483 

2b)  }•!„«*  1.0,  PROPAGATING  MODE 


a 

c 

r 

Ci 

,10 

.824508 

.095961 

15 

.790995 

.094567 

,20 

.766441 

.090119 

,25 

.746227 

.083887 

,30 

.728650 

.075980 

,35 

.713325 

.066800 

,40 

.700178 

.056500 

45 

.689115 

.045423 

.50 

.680052 

.034040 

60 

65 


667009 

662642 


011367 

000350 


a 

.10 

.15 

.20 


2c)  m  =1.2,  PROPAGATING  MODE 


.808530 

.768575 

.746423 


.25 

.30 

.35 


.728869 

.713971 

.701157 


.40 

.45 

.47 


.6^9315 

.681388 

.678352 


.007919 

.067700 

.060019 

.051536 

.041967 

.031200 

.019372 

.006778 

.001596 


2d)  M^O,  STANDING  MODE 


a 

c 

r 

ci 

.1 

0.0 

.231871 

.2 

0.0 

.351572 

.3 

0.0 

.421749 

.4 

0.0 

.467320 

.5 

0.0 

.495630 

.6 

0.0 

.509976 

.7 

0.0 

.512729 

.8 

0.0 

.506042 

.9 

0.0 

.491820 

1.0 

0.0 

.471650 

1.1 

0.0 

.446790 

1.2 

0.0 

.418208 

i  •  3 

0.0 

. 386635 

1.4 

0.0 

.352607 

1.5 

0.0 

.316510 

1.6 

0.0 

.278608 

1.7 

0.0 

.239067 

1.8 

0.0 

.197972 

1.9 

2.0 

2.1 


0.0 

0.0 

0.0 


155345 

111149 

065246 


2  4 


.70 

.80 

.90 

1.00 

1.10 

1.20 

1.23 


2e)  M<n=  1.0,  STANDING  MODE 
c 


0.0 

.325303 

0.0 

.376416 

0.0 

.434305 

0.0 

.454416 

0.0 

.450267 

0.0 

.429563 

0.0 

.396980 

0.0 

.355153 

0.0 

.305120 

0.0 

.246296 

0.0 

.175520 

0.0 

.081841 

0.0 

.043534 

O'  > 


2f )  H^-1.4,  STANDING  MODE 


1 


.1  0.0  .26/995 
.2  0.0  .395734 
.3  0.0  .441200 

0.0  .441578 
0.0  .414618 
0.0  .368707 


7 

8 
9 


0.0 

0.0 

0.0 


.306499 

.224023 

.093603 


2g)  11^=2. 0,  STANDING  MODE 


a 


Ci 


.15 

.20 

.30 


0.0 

.381035 

0.0 

.  A25659 

0.0 

.444A03 

.AO 

.50 

.60 


0.0 

.409147 

0.0 

.339340 

0.0 

.229747 

• 

65 


0.0 


137767 


00  n|  Ot  V' 


INSTABILITY 


CHARACTERISTICS  FOR  THE  ASYMMETRIC  JET 


TABLE  3. 


a 


.1 

.2 

.3 


1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

2.5 

3.0 

3.5 
4.0 


3a)  M  =0.0,  SINUOUS  MODE 
*  00 


c 

r 

Ci 

.022716 

.057622 

.096191 

.034685 

.075168 

.106869 

.134367 

.169211 

.199485 

.126540 

.135130 

.134785 

.223771 

.241919 

.254235 

.127902 

.116944 

.104212 

.261494 

.265033 

.266183 

.091611 

.080302 

.070682 

.265879 

.264760 

.263194 

- 

.062672 

.056016 

.050437 

.261399 

.259499 

.257556 

.045706 

.041642 

.038108 

.255643 

.253755 

.245207 

.035004 

.032250 

.022082 

.238341 

.233025 

.229011 

.015553 

.011077 

.007910 

I 


VO  OS  >4  O'  Ul 


3b)  1.0,  SINUOUS  MODS 

c 


2  S 


.i 

.2 

.3 


1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

2.5 

3.0 

3.5 
4.0 


.025250 

.066401 

.114810 

.163526 

.206953 

.240919 

.263146 

.271909 

.271495 

.268354 

.264800 

.261773 

.259106 

.256705 

.254633 

.252696 

.250912 

.249280 

.247718 

.246241 

.239841 

.234696 

.230608 

.227434 


.038126 

.083135 

.116629 

.132723 

.132538 

.119500 

.098537 

.076103 

.059204 

.048282 

.041100 

.036078 

.032331 

.029392 

.026992 

.024966 

.023213 

.021664 

.020275 

.019014 

.014013 

.010404 

.007697 

.005654 


VO  OO  "-l  O'  c/i  o- 


3c)  M  -2.0,  SINUOUS  MODE 

U1 


1 


c 

Ci 

a 

r 

2 

.092437 

.102199 

3 

.169565 

.130870 

4 

.242842 

• 

.123709 

5 

.315747 

.081390 

,6 

.405591 

.064530 

,7 

.470175 

.050658 

.8 

.525611 

.040951 

,9 

.571746 

.036654 

0 

.607611 

.035926 

1 

.633812 

.035146 

,2 

.652273 

.035042 

.3 

.665219 

.031993 

.4 

.674078 

.027458 

,5 

.679530 

.022257 

.6 

.682270 

.017361 

,7 

.683250 

.013417 

.8 

.683357 

.010509 

.9 

.683125 

.008409 

,0 

.682793 

.006866 

,5 

.681431 

.003002 

3d)  Mw=3.0,  SINUOUS  MODE 


a 

c 

r 

ci 

.1 

.041923 

.058484 

.2 

.134351 

. 120942 

.3 

.250400 

.125258 

.4 

.364384 

.112313 

.5 

.441116 

.093284 

.6 

.508862 

.070646 

.7 

.568281 

.057581 

.8 

.613491 

.053023 

.9 

.645471 

.049900 

1.0 

.668287 

.045140 

1.1 

.684931 

.037703 

1.2 

.696740 

.026809 

1.3 

.700800 

.010626 

1.4 

.692436 

.002767 

1.5 

.688292 

.001691 

1.6 

.686144 

.001325 

1.7 

.684820 

.001130 

1.8 

.683920 

.001000 

1.9 

.683265 

.000899 

2.0 

.682766 

.000814 

3e)  M  =4.0,  SINUOUS  MODE 

to 


.1 

.054213 

.071056 

.2 

.185692 

: 129058 

.3 

.343024 

.122561 

.4 

.440501 

.105400 

.5 

.523520 

.078015 

.6 

.590604 

.067213 

.7 

.636500 

.061840 

.8 

.668800 

.055668 

.9 

.692983 

.046829 

1.0 

.713555 

.034536 

1.1 

.736208 

.020525 

1.2 

.760259 

.011891 

1.3 

.780733 

.007100 

1.4 

.798010 

.003777 

1.5 

.812881 

.001170 

^  m  vO  IN  flO  ON 


.1 

.2 

.3 


1.0 

1.1 

1.2 


1.3 

1.4 


1.5 


1.6 

1.7 

1.8 


1.9 

2.0 

2.5 


3.0 

3.5 

4.0 


3f)  M  =0,  VARICOSE  MODE 
00 

c 


.791496 

.694068 

.634863 


.596700 

.572033 

.556473 


.548758 

.548309 

.553977 


.564245 

.577329 

.591419 


.605205 

.617916 

.629175 


.638852 

.646961 

.653605 


.658945 

.663174 

.673978 


.677488 

.678856 

.679452 


176117 

187158 

,173131 


,152443 

,130738 

,110405 


.092825 

.078848 

.068563 


.061258 

.055838 

.051351 


.047201 

.043095 

.038880 


.034763 

.030638 

.026674 


.022975 

.019621 

.008493 


.003768 

.001753 

.000846 


3^ 


3g)  M  -1.0,  VARICOSE  MODE 


a 


3h)  M  =1.3,  VARICOSE  MODE 

CO 

c 

r 


3  + 


c 


i 


.01 

.05 

.10 

.15 

.20 

.30 


. 888249 
.793824 
.  727671 

.685409 

.656887 

.628048 


045163 

096293 

103920 

095711 

080776 

041152 


35 

40 


.628602 

.634545 


016927 

005976 


9M 


Figure  20b 


Figure  21a  Moo  "  1.0  a  =  .90  c  =  .480444+. 1615691 


Figure  21b 


Figure  23a  M®  ■  1.0  a  m  .30  c  ■  . 728650+. 0759801 


Figure  23b 


Figure  27b 


